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Porcupines are networks of gravitational wave detectors in which the detectors and the distances 
between them are short relative to the gravitational wavelengths of interest. Perfect porcupines 
are special configurations whose sensitivity to a gravitational plane wave is independent of the 
propagation direction or polarization of the wave. I develop the theory of porcupines, including 
the optimal estimator h 13 for the gravitational wave field; useful formulae for the spin-averaged and 
rotationally-averaged SNR 2 ; and a simple derivation of the properties of perfect porcupines. I apply 
these results to the interesting class of "simple" porcupines, and mention some open problems. 
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This decade, we hope and expect that gravitational 
waves will be directly detected for the first time; this will 
mark the beginning of what promises to be a long and 
fruitful era of gravitational wave astronomy. One can 
only guess what wonders will be revealed when this new 
window onto the universe is flung fully open. As atten- 
tion shifts from gravitational wave detection to gravita- 
tional wave astronomy, so will it shift from individual 
gravitational wave detectors to networks of multiple de- 
tectors that that function together as gravitational wave 
telescopes. (For previous work on gravitational wave net- 
works, see [l|-|4| and references therein.) There is an im- 
portant regime in which the individual detectors and the 
distances between them are short relative to the gravi- 
tational wavelengths of interest, so that the arms of the 
various detectors in the network may be thought of as 
emanating from nearly the same point in space, like the 
fanned quills of a frightened porcupine. In a previous pa- 
per [U , such "porcupines" were motivated by recent ideas 
and developments in gravitational wave detection 
but it is worth reframing the motivation in more gen- 
eral terms: one may suspect that porcupines will play 
an important role in gravitational wave astronomy for 
the simple reason that astrophysical gravitational wave- 
lengths are typically very long. One will often be inter- 
ested in gravitational wavelengths that are longer than 
the practical size of one's network (e.g. the diameter of 
the Earth, or the Solar System). 

Ref. [1| focused on "perfect porcupines": special con- 
figurations with the defining property that their sensi- 
tivity to an incident gravitational plane wave is inde- 
pendent of the propagation direction or polarization of 
the wave. (A number of other nice properties follow 
automatically once this condition is satisfied.) Here I 
develop a more natural and general formalism for han- 
dling porcupines (perfect or imperfect). This formal- 
ism yields, for an arbitrary imperfect porcupine: (i) the 
minimum- variance unbiased estimator h %3 for the gravita- 
tional wave field; and (ii) useful expressions for the spin- 
averaged and rotationally-averaged SNR 2 . It also yields a 
simple derivation of the properties of perfect porcupines. 
I apply these results to the interesting class of "simple" 
porcupines (defined below) and discuss open problems. 



GENERAL FORMALISM 

To fix notation, my fourier conventions are 

/•oo poo 

9(t)= dfg(f) C -^f\ ~g(f)= dtg(t)e +2 ^ t . (1) 



The lower case latin indices {i, j, k, I, m, n} label the 3 
spatial directions: i, j, k, l,m,n = 1, 2, 3. The upper case 
latin indices {A, B} label the 2 gravitational wave polar- 
izations: A, B = 1,2. The lower case greek indices {a, /?} 
label the N detectors in the network: a, f3 = 1, . . . ,N. 
Indices in square braces [. . .] are symmetrized. I use the 
Einstein summation convention: repeated indices (one 
upper, one lower) are summed. I use hats in two differ- 
ent ways: (i) to denote variables with a random noise 
contribution (such as estimators); and (ii) to denote unit 
3- vectors. I have tried to make the context clear enough 
to avoid confusion of these two meanings. 

A gravitational wave on Minkowski space is described 
in transverse-traceless gauge [TlJ] by the line element 



ds 2 



-dt 2 + [Sij + 2h^\dx % dx j . 



(2) 



In this metric, worldlines with x = constant are 
geodesies; along such worldlines, the proper time t co- 
incides with the coordinate time t. When gravitational 
waves reach us from a distant astronomical source, they 
appear to us as plane waves travelling in the n direction 



h«(\)=P%(n)h A (\) 



(3) 



where the polarization waveforms h A (X) are arbitrary 
functions of A = t — n ■ x, and the polarization tensors 
P% (n) form an orthonormal basis on the 2-dimensional 
space of symmetric, transverse, traceless 3x3 matrices: 



hiPjl (n) 

pmhyp^h) 



0, 
0, 
0, 

6 B , 



(4a) 
(4b) 
(4c) 
(4d) 



Let us imagine that all of the detectors are situated at 
x = 0. The estimator h a represents the measured output 
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of detector a, a sum of genuine gravitational wave signal 
h a (t) and noise n a {t): 



h a (t) = h a (t)+n a {t). 



(5) 



If we assume the network's response is linear and time- 
translation invariant, we can write 



h a (t) 



+ 00 



dTW t a JT)h lj (t-T). 



(6) 



Since h a (t) is real and h l ^(t) is real and symmetric, 
W£j(t) is also real and symmetric. We model the noise 
n a (t) as stationary and gaussian, with zero mean, so it 
is characterized by its correlation function C a ^(T) or, 

equivalently, its spectral density S a /3 (f) = C a o(f): 



(h a (t + T)n p (t)), (7a) 
(n a (/)W)>- ( 7b ) 



S a o{f) induces a natural inner product on the space of 
signals (or noise) in the network: 

/oo 
<VPM)[s-Hf)] a f>i p (f)- (8) 
-00 

Using matched filtering, a given gravitational wave signal 
h % i(t) can be detected by such a network with expected 
signal-to-noise (SNR) given by 



SNR 2 = (ft|ft) : 



tfhMlS-HfWpWif) (9a) 

O 

dfh* 3 (f)K^ kl (f)h kl (f) (9b) 



where we have introduced a kernel 



K ij M (f) w^ms-HfTp K(f) (10) 



with the following properties 



K 



-I j 



uU)=K [i \ kl] (f)=K^{fr=K^ kl {-fr 



(ii) 



If a gravitational wave signal (which depends on various 
parameters fi a ) is detected, and the likelihood function 
may be approximated as gaussian cx exp[— (l/2)/^ a r a (,/j b ] 
near its peak, then the expected inverse covariance ma- 
trix is the Fisher information matrix, given by 



r„ 



/ dh dh 



(12) 



Now consider the angular resolution of this network. Let 
us define an orthonormal triad consisting of n and two 
additional unit vectors rhp (fi = 1,2). Under a rotation 
around the direction rhp by an infinitcssimal angle 9^ , 
h ij (t) transforms as -> R\R\h kl (t) where R l 3 

6 % j — e i j k rhp0' i . So the angular part of F a b is given by 



/ dh 



'Hi 



where we have defined 



(14a) 
(14b) 



For fixed /2 and D, L l ?- kl (f) has properties exactly akin 
to those in (fTT|) . The "off-diagonal" terms in the Fisher 
matrix, with one index p, corresponding to an angular 
parameter 9^, and the other index \ corresponding to 
any other (non-angular) parameter is 



50 Flh* ( f\ 

df^^P^hYK l \ l {f)e\ s rnlh^{f). (15) 

00 



OPTIMAL ESTIMATOR FOR h %1 

We want to find the minimum-variance unbiased esti- 
mator ft SJ '(/) and its covariance. Since h a (f) is linearly 
related to [h a (f) = (/)/#(/)], ft«(/) is lin- 

early related to h a {f): 

W(f) = V i i<J)h a (J). (16) 
It is stationary and gaussian, with covariance 

(Sh ij (f)6hUf)) = S ij kl (f)S(f - f) (17) 

where Sh ^ = h 3 — h % K Our task is to express V 1 ^ and 
S lJ ' kl in terms of the given tensors W 7 ," and S a g that define 
the porcupine. 

The estimator is unbiased if (ft y ) = ft*- 7 ', and hence 
h ij = V%W%h kl . Since this must hold for arbitrary h ij , 
we learn that 



where 



kl 



6 [i S i] 
[k° i] 



(18) 



(19) 



is the identity operator on the 6-dimcnsional vector space 
of 3 x 3 symmetric matrices. The estimator hij only exists 
if the porcupine tensor W£- has a left-inverse V l J in the 
sense of (|18J). In a porcupine consisting of N detectors: 
if N < 6, V % i docs not exist; if N = 6 then, when V% 
exists, it is specified uniquely by (JT5J; and if N > 6 then, 
when V Z J exists, it is not specified uniquely by ((T8j) . In 
this last case, the residual ambiguity in V 1 ^ may be fixed 
by using the method of Lagrange multipliers to minimize 
the total variance P ] kl S \j of the estimator fty, subject 
to the constraint (fl8|) : for help, see e.g. Appendix D in 
[l3T ] . In this way we find 



v i 2 = (D- 1 ) ii M wt(s- 1 )f > a 



and 



kl 



(20) 



(21) 
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where we have defined the tensor D and its inverse D 1 



D' 



kl 



wris-^w^, 



l\mn 



kl 



(22a) 
(22b) 



This completes our derivation of the minimum- variance 
unbiased estimator K lJ and its covariance S lJ k ,. Let us 
add a few remarks. 

Remark 1. The minimum- variance unbiased esti- 
mator for the trace of h lJ is h = Sijh*-' , with variance 

(kf)km - (HiMkm = h^uU^u - n 

If the porcupine detects a gravitational wave, h should be 
consistent with zero, to within this predicted uncertainty; 
this is an important observational test to distinguish gen- 
uine gravitational waves from spurious signals. 

Remark 2. If 8 %3 W{j = (as in a network of equal- 
arm Michelson interferometers like LIGO/ VIRGO 01) 
the porcupine is insensitive to the trace of h lJ , and we 
can only construct an estimator for the traceless part of 
h' lJ . It is convenient to introduce the notation [. . .] T to 
denote the traceless part of a tensor. We must be careful, 
since this paper involves the relationship between several 
different linear spaces, each with its own trace. In par- 
ticular, [. . ] T will denote tracelessness with respect to 
3-dimensional traces {i.e. those that can be performed 
using Sij and 5 lJ , but not <5 4 -). Explicitly: 



[T ij ] T =T ij --S ij S k lT kl 



(23a) 



[T 1 



k-il 



1 



kl 



W kl 



Tg« Vpq-L mn U ffc 



(23b) 



Thus, we want the minimum- variance unbiased estimator 
[h i3 ] T and its covariance. This is done by following the 
same steps as above; in fact, the answer is still given 
by Eqs. ([16] -[22]), as long as we replace each factor in 
these equations by its traceless part: e.g. h lJ — > [/i lJ ] T ', 



ki 



kit 



Here 



tjij IT 
I 1 Ml 



s [i s j] 

[k° I] 



~-8 lJ 5 



(24) 



is the identity operator on the 5-dimensional space of 
3x3 symmetric traceless matrices. The estimator [h lJ ] T 
only exists if W°- has a left- inverse : } T in the sense of 
(the traceless part of) Eq. (fT8|) : and this is only possible 
for TV > 5 detectors. 

Remark 3. Although the porcupine may be simul- 
taneously immersed in many gravitational plane waves, 
travelling in many different directions, suppose that there 
is at most one plane wave in each frequency bin f ± 5f, 
traveling in an unknown direction n(f). Then, since the 
gravitational wave is transverse, we have h lJ (f)nj(f) 
!). which implies that Det[h lJ (f)] should vanish. We can 
check this observationally by confirming that the esti- 
mator Det[/i lJ (/)] is consistent with zero, to within the 



predicted uncertainty. (Calculating the predicted vari- 
ance of this estimator is an exercise in applying Wick's 
Theorem.) Also note that if we choose two mutually 
orthogonal unit vectors a and b, and use these to form 
the estimators A 1 = h l3 a,j and B l = ti lJ bj, then N = 
(A x B)/\ A x B\ is an estimator for the zero-eigenvector 
of hij {i.e. ±n, the propagation direction up to a sign). 



AVERAGE SNR 2 

Let us average the expected SNR 2 ©: first over the 
polarization state, and then also over the direction h of 
the incident gravitational wave. First we average over 
the polarization, keeping h fixed, to obtain: 



where 



SNR? 



dflkf^v^m^kiif) 

) 

ft«(/)ft* J (/), 

(l/2)jVlW fc 1(n)* 
id 



(25) 

(26a) 
(26b) 



= (l/^A^A^, - (l/4)A«A Wj (26c) 
and Ay- = <S y - — fiifij is the projector onto the plane 



perpendicular to h. Now use Eq. (2.3b) in [12| to find 
J^ ij k i(n) = l[I ij kl f (27) 

and use this identity to obtain the SNR 2 , fully averaged 
over both the polarization and direction n of the the in- 
cident gravitational wave: 

" 1' 



(SNR 2 



df\h{f)\--[I kl l3 ] T K^ kl {f). 



(28) 



PERFECT PORCUPINES 

A perfect porcupine is a network for which the tensor 
K 1J kl {f) is rotationally invariant. The only objects from 
which we can build such a tensor are 5 iJ and e tlk . Im- 
posing the constraints (fTTI) we find that the most general 
perfect porcupine is of the form 

K ij ki(f) = FU)[I l \i\ T + G{f)8 l3 5 kl (29) 

where F and G are real and symmetric under / — > — / 

F{f)=F{f)* = F(-f), G{f) = G{f)* = G(-f). (30) 

Then the expressions for the SNR 2 , and the angular parts 
of the Fisher matrix, simplify beautifully 

SNR 2 =/_ 00 co d/F(/)|M/)| 2 , (31a) 

i> = (snr 2 )^ p , r px = r x ^ = o. (3ib) 

Thus, a perfect porcupine localizes a gravitational wave 
source within a circular spot (really two antipodal spots) 
of radius 1/SNR, regardless of the source's properties. 
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SIMPLE PORCUPINES 

Suppose that the individual detectors in the network 
are identical to one another (up to spatial orientation), 
with noise that is uncorrelated between different detec- 
tors: S a j3 (f) = S(f)S a p. Further suppose that each de- 
tector's frequency response may be factored out from its 
tensor structure: W$(f) = W(f)A?j. (Wc can take 
the frequency-independent matrices Afj to be normal- 
ized: AQ* Afj = 1 with a unsummcd.) A network with 
these properties will be called "simple." 

The averaged SNR 2 (|28p for a simple porcupine is 



(SNR 2 



N 



(3- \TrA\ ) 
15 



mimnl (32) 

, 1 s(f) (32) 



PA*. 



where Tr A 

A "simply perfect porcupine" is both simple and per- 
fect. In this case, we can trace Eq. (|29p in two incquiv- 
alent ways (by contracting with 6 ' \S j or o~ij5 kl ) to find 



F(f) = N 
G(f) = N 



\W(f)\ 2 (3-\TvA\ z ) 

S(f) 
\W(f)\ 2 \TrA 

S(f) 9 



15 

i2 



(33a) 
(33b) 



For such a porcupine, we see that the SNR 2 expression 
(|31a|) coincides with the averaged expression (|32|) for a 
simple but imperfect porcupine. Note that these expres- 
sion depend only on N, S(f), W(f), and Tr A; and not 
on the directional orientation of the individual detectors 
that make up the network. In other words, if two dif- 
ferent simply perfect porcupines are built from the same 
collection (i.e. the same number and type) of individual 
detectors, then they will have identical properties, de- 
spite the fact that they may seem like entirely different 
network configurations. 

Now consider the problem of finding simply perfect 
porcupine configurations explicitly. Two cases deserve 
particular attention. In the first case, the network is 
constructed from single- arm detectors (like the proposed 
AGIS detector 0S,L3LQl); the arm of detector a lies along 
the direction k a so that A% = kfkf and TrA = 1. 
In the second case, the network is constructed from 
double-arm detectors (like the LIGO/VIRGO Michel- 
son interferometers 0); the two arms of detector a lie 
along the two perpendicular directions p a and q a so that 
Ag. = (pfpj - q?qf)/\/2 and TrA = 0. In these two 
cases, Eq. ([29]) becomes 



T, a kiki^k l a = (N/15)(d»S kl +S ik S^ l +S u S^), (34a) 
£„ ~ ) {PtPf ~ « ) = (2JV/5) [I ij kl ] T • (34b) 



{i,j, k, I}; this gets us down to 15 independent equations, 
one of which (obtained by tracing both sides with 5ij5u) 
is automatically satisfied. Thus we have 14 independent 
equation for 2N — 3 variables (2 angles for each k a , minus 
3 angles corresponding to an arbitrary rigid rotation of 
the entire network). So, for N > 9, Eq. (|34a[) will have a 
(2N — 17) parameter family of incquivalent solutions. In 
addition, there is a unique smaller solution with N = 6 
detectors along the 6 diameters of an icosahedron [1, [H [f| . 
On the other hand, for the double-arm Eq. (|34b[) . there 
are no solutions for N < A, and continuous families of 
incquivalent solutions for all N > 5. 



DISCUSSION 

Let us highlight three open problems, (i) The first 
problem is to find the general solution of Eq. (|34ap or 
Eq. (|34b[) [that is, the most general configuration of a 
(single-arm or double-arm) simply perfect porcupine]. 
Which among these configurations could be most prac- 
tically situation on the Earth's available land area? (ii) 
Whereas perfect porcupines have a completely isotropic 
sensitivity pattern, "pointed porcupines" lie at the oppo- 
site extreme: they are the porcupines with the sensitiv- 
ity pattern that is most highly localized on the sky. The 
second problem would be to develop the theory of such 
pointed porcupines, particularly since they may be the 
configurations that maximize the expected event rate (for 
a fixed collection of detectors, and assuming the gravi- 
tational wave sources are Poisson distributed through- 
out space). The third problem is to extend the present 
formalism in a relativistically correct way to investigate 
networks in which each detector is allowed to follow its 
own (possibly accelerating and rotating) trajectory. 
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Note that both sides of the single-arm Eq. (|34ap arc 
completely symmetric under permutation of the indices 
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